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The quantum chaos in the finite-temperature Yang-Mills-Higgs system is studied. The energy 
spectrum of a spatially homogeneous SU(2) Yang-Mills-Higgs is calculated within thermofield dy- 
namics. Level statistics of the spectra is studied by plotting nearest-level spacing distribution 
histograms. It is found that finite temperature effects lead to a strengthening of chaotic effects, 
i.e. spectrum which has Poissonian distribution at zero temperature has Gaussian distribution at 
finite-temperature. 
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I. INTRODUCTION 

Quantum chaos is a relatively new area in physics and 
has been the subject of extensive studies for the last two 
decades It has found applications in atomic and 

molecular physics, nuclear physics and condensed matter 
physics. In the past few years there is growing interest in 
quantum chaos in particle physics, too. Being the quan- 
tum theory of classically chaotic systems the quantum 
chaology studies fluctuations in the energy spectra and 
wave functions of such systems. It is well known, that the 
energy spectra of systems, whose classical counterparts 
are chaotic, has the same statistical properties as those 
for random matrices. Therefore one of the main topics 
in the quantum chaology is to study the statistical prop- 
erties of the classically chaotic systems. Recently energy 
fluctuations and quantum chaos in hadrons and QCD has 
become a subject of theoretical studies In partic- 

ular, it is found that the quark-gluon system in QCD is 
governed by quantum chaos in both confined and decon- 
fined phases [H, Q ■ The statistical analysis of the mea- 
sured meson and baryon spectra shows that there is quan- 
tum chaos phenomenon in these systems . The study 
of the charmonium spectral statistics and its dependence 
on color screening has established quantum chaotic be- 
haviour 7}. It was claimed that such a behaviour could 
be the reason for suppression 0. 

In recent years there has been considerable interest to 
determine the role of dynamical chaos in field theories 
0-123 . Chaoticproperties of Yang-MillsEl , Yang- 
Mills-Higgs an d Abelian Higgs |[20| have been 
treated. The main point in this considerations is the fact 
that the Hamiltonians of the Yang-Mills and Yang-Mills- 
Higgs systems can be written in the same form as those 
for the coupled nonlinear oscillators. This allows one 
to use, for their treatment the same methods as in the 



case of coupled nonlinear oscillators. Quantum chaos in 
Yang-Mills-Higgs system was also studied recently [2^, 
|23| . However all works on chaos in field theories and 
hadrons are restricted to considering zero-temperature 
cases. 

In this paper we study quantum chaos at finite- 
temperature in the Yang-Mills-Higgs system. Recent ad- 
vances in heavy ion collision experiments allow one to 
create hot and dense hadronic and quark-gluon matters. 
The role of finite-temperature effects in such systems be- 
comes important. Especially in quark-gluon or nuclear 
matter hadrons behave as complex systems where strong 
level fluctuations can be observed. Thus the role of finite- 
temperature effects as well as level fluctuations are cru- 
cial in such a systems, that obviously leads to a need for 
studying quantum chaos at finite-temperature. In partic- 
ular, finite-temperature effects cause fluctuations in their 
energy spectra. 

We address the problem of heat-bath in quantum chaos 
through the thermofield dynamics (TFD) formalism, a 
real time finite temperature field theory [2J, |2|| |2(| |23, 
|2S| . TFD is interesting for our proposal here by its 
remarkable algebraic structure (this is not the case of 
the Matsubara [29j or the Schwinger-Keldysh [3(| for- 
malisms). Actually the central ideas of TFD involve an 
algebraic doubling in the degrees of freedom and a Bogoli- 
ubov transformation giving rise to the thermal variables. 
As we will show, this TFD prescription is a useful tool 
to explore finite-temperature effects in the energy fluctu- 
ations of the Yang-Mills-Higgs Hamiltonian, represented 
in terms of the annihilation and creation operators. 

In section II we present the Yang-Mills-Higgs system 
at zero temperature. In section III the Yang-Mills- 
Higgs system at finite-temperature is studied using the 
TFD formalism. In section IV we present numerical re- 
sults for finite-temperature and compare them to zero- 
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temperature results, thus bringing out crucial role of 
finite-temperature in Yang-Mills-Higgs system and quan- 
tum chaos. Finally in Section V we present some conclu- 
sions and new directions to pursue the important ques- 
tion of quantum chaos. 

II. ZERO TEMPERATURE CASE 

The Lagrangian for Yang-Mills-Higgs system with 
SU(2) symmetry is given as 

L = -\f« v F^ + 1 -{D,4>) + (D^) - V(<f>) (1) 

where 

F^ = d„A a v -d v A^+gAlA c v 

(D^) = drf - igA\T h <t> 

with T b = a/2, 6=1,2,3 generators of the SU(2) alge- 
bra and g is a coupling constant. The potential of the 
scalar (Higgs) field is 

F(0)^ 2 H 2 + AH 4 , 

where p, and A are constants. Here we give a brief descrip- 
tion of the non-thermal case [2^|. In (2 + l)-dimensional 
Minkowski space and for spatially homogeneous Yang- 
Mills and Higgs fields which satisfy the conditions 

5^ = 0^ = 0, i =1,2,; 

and in the gauge Aq = 0, the Lagrangian can be written 
as 

L = $ 2 + \(A\ + A 2 ) - g[\A 2 A 2 \(A 1 A 2 ) 2 + 

(A 2 + - {Ai4>) 2 ~ [A 2 $f] - V($), (2) 

where $ = {<p\ c/) 2 , 3 ), A x = (A\,A\,A\) and A 2 = 
(A\, A 2 ,, A 2 ) 

To treat the chaotic dynamics it is convenient to use 
the canonical formalism and work with the Hamiltonian 
of the system instead of the Lagran gian . The Hamilto- 
nian of the system can be written as 22] 

H = \(p\ + Pi) + 9 2 v 2 (ql + <£) + \ 9 V iq l (3) 

where <fio = (0, 0, v) q% = A\, q 2 = A 2 , (other components 
of the Yang-Mills fields are zero) p\ = qi and p 2 = q 2 , 
with cu 2 — 2g 2 v 2 being the mass term of the Yang-Mills 
field. This is the Hamiltonian of the classical system. 
Replacing pi and g, by operators and introducing the 
following creation and destruction operators 




we obtain the corresponding quantum Hamiltonian: 

H = H + X -g 2 V (4) 

where 

H = Lu(aia~l + a 2 a 2 + 1), 

and 

V = ^ 2 -(a 1 +a+) 2 (a 2 + a+) 2 , 

with tu 2 = 2g 2 v 2 and the operators du and a~l satisfy 
the commutation relations [cife,^] = Ski, k, I = 1,2. 
The eigenvalues of this Hamitonian are calculated by nu- 
merical diagonalization of the truncated matrix of the 
quantum Yang-Mills-Higgs Hamiltonian in the basis of 
the harmonic oscillator wave functions p2| . The matrix 
elements of Hq and V are 

< n' 1 ,n 2 \H \ni,n 2 >= w(ni + n 2 + l)5 n / ni 5 n / na) 

and 

< n' l ,n' 2 \V\nx,n 2 >= -^{y/n^m - l)^n' 1 m-2+ 

\/m(m - l)S n ' ini+2 + (2m + l)6 n > ini } 
*{\/n 2 (n 2 - l)<5„^ n2 _2 + 



\/n 2 (n 2 - l)<5„ 2 n 2+2 + (2n 2 + l)S n > 2 „ 2 } 

The numerical procedure for diagonalization of this ma- 
trix is described by Salasnich [23]. We use the same 
method in the case of finite-temperature calculations. 

III. FINITE-TEMPERATURE CASE 

To treat quantum chaos in the finite-temperature 
Yang-Mills-Higgs system, we apply Thermofield Dynam- 
ics(TFD). TFD is a real time operator formalism of quan- 
tum field theory at finite temperature in which any physi- 
cal system can be constructed from a temperature depen- 
dent vacuum |0(/3) > which is a pure state. The thermal 
average of any operator is equal to the expectation value 
between the pure vacuum state |0(/3) > defined by ap- 
plying Bogolyubov transformations to the usual vacuum 
state . Furthermore, TFD has two main features. The 
first one is the doubling of the Fock space such that the 
original Fock space and its double are defined non-tilde 
and tilde space respectively. All operators are also dou- 
bled and the finite-temperature creation and annihilation 
operators are constructed by Bogoluybov transformation 
between tilde and non-tilde operators. This is the same 



3 



procedure in writing down the vacuum state at finite tem- 
perature. Mathematically, the field operators have the 
following properties: 

(AiAj) — AiAj , 
(cA i + A j ) = c*A i + A j , 

(4*) = (4)+ 
(A) = Ai 

[MAj] = 0. 

The Yang-Mills-Higgs Hamiltonian in TFD is given as 
H = H-H (5) 
where H is given by eq0] and H is given as 

H = H Q + X -g 2 V 

with 

H = u)(aia^ + fi 2 fij + 1), 

and 

V=^(ai + ~a+) 2 (~a 2 + a+) 2 , (6) 

It has been established [3lJ that in an algebraic ap- 
proach the doubled set of operators may be considered as 
a set of operators that relate to the physical observables, 
O, and a second set that are generators of symmetries, O. 
The hat operators are responsible, in particular, for time 
development and are needed for the purpose of scatter- 
ing, decay and any transitions between states. The phys- 
ical observables lead to the quantities that are measured 
in experiment. Both for physical observables and gen- 
erators of symmetry after Bogolyubov transformations 
leading to finite-temperature creation and annihilation 
operators and to a pure vacuum state only the non-tilde 
operators are required for getting the appropriate matrix 
elements. However, it is clear that for an analysis of any 
system at finite temperature both set of operators , O 
and O, are needed since it is important to generate the 
appropriate symmetry of the system in any time develop- 
ment while considering the matrix element we have de- 
cided to investigate the chaotic behaviour for the Hamil- 
tonian that relates to physical observables, H and the 
Hamiltonian that generates the symmetry of the system, 
H. Results of the finite-temperature quantum chaos in 
Yang-Mills-Higgs system will be displayed for the physi- 
cal observables, and for the generators of symmetry, H. 
Only then we will draw conclusions about the approach 



of quantum chaotic behaviour in the finite-temperature 
quantum field theory case of Yang-Mills-Higgs theory. 

First wc need to rewrite the non-tilde part of the 
Hamiltonian in the temperature-dependent form using 
the Bogolyubov transformations which are given by 

flfc = a k ((3)cosh0 + a^(j3)sinh9 
a t = a~^{(3)cosh9 + a k (j3)sinh9 
fife = ((3)sinh6 + a k ((3)cosh6 
at = a k ((3)sinh9 + at (/3)cosh6 

where 



where tilde and non-tilde creation and annihilation oper- 
ators satisfy the following commutation relations: 

[a k ((3),a+(f3)} = S kl [a k ((3),a+(0)] = 5 U 

l,k = 1, 2, and sinh 2 9 = (e 13 — l) -1 . All other commuta- 
tion relations are zero. 

Then the temperature-dependent forms of Ho and Ho 
are 

H =uj{(F 1 +F 2 )cosh 2 e+ 

(Li + L 2 )sinh 2 8 + (Si + SzjcoshB 'sinh9 + 1}, 

H = uj{{F 1 + F 2 )sinh 2 9+ 

(L 1 + L 2 )cosh 2 e + (Si + S 2 )coshOsinhQ + 1}, 
where 

L k = a k (0)a+({3), 

S k =a+((3)~a+(l3) + a+((3)a k (P), 
For V and V we have 

V = -^{AxcostfO + B 1 cosh0sinh6+ 

Cisinh 2 0}{A 2 cosh 2 9 + B 2 cosh9sinh8 + C 2 sinh 2 0} 

V = -J-^{A 1 sinh 2 e + B x cosh0sinh6 + Cicosh 2 6}x 
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x{A 2 sinh 2 9 + B 2 cosh9sinh9 + C 2 cosh 2 9}, 

where 

A k = {a k {p)+a+{p))\ 

B k = {a k {f3) + a+ (/?))(£+ (/?) + a k (/?))+ 
(a+(/3) + a fe (/3))(a fe (/3) + a+(/J)), 

and 

C fc = (a+(f3)+a k (P)) 2 . 

In the first approach the energy eigenvalues of the ther- 
mal Yang-Mills-Higgs system can be calculated by di- 
agonolization of the following matrix: 

R =< n[n' 2 , n[ h' 2 \ H + ^g 2 V\n 1 n 2 ,n 1 h 2 > . (7) 
The elements of the matrix can be calculated explicitly: 

< n'-^n^, n' 1 n 2 \HQ\nin 2l h\h 2 >= 

= uj{{ni +n 2 + 1)(1 + 2sinh 2 9)6 n > ini 5 n > 2n2 + 

+(ni+l)S n ' i7ll+ iS n ' 2n2 + 
(n 2 + l)d n ' 2 n 2 +iS n ' ini )}cosh9sinh9} 

and for V 

<ri 1 ri 2l f)! 1 f)! 2 \y\n\n 2 ^h\h 2 > = 

— ^ < n[n' 2 ,h' 1 h' 2 \{A 1 A 2 cosh 4 9 + dC 2 sinh 4 9+ 

(AiC* 2 + BiB 2 + C 1 A 2 )cosh 2 9sinh 2 9+ 

{A X B 2 + B 1 A 2 )cosh 3 9sinh9+ 

{B 1 C 2 + CiB 2 )cosh9sinh 3 9}\nin 2 ,hi,h 2 > (8) 
The matrix elements of A k , B k and C k are given as 

< n[n' 2 ,h[h' 2 \A k \nin 2 ,nih 2 > = 

= ^n k (n k - l)S n > k , lk _ 2 + (2n k + l)<5 n ^„ fc + 



< n\n 2 ,h[n' 2 \B k \nin 2 ,nih 2 > 
= 2nk<>n k n k -i + (2n k + l)S n > kUk+1 , 

and 

< n' l n' 2 ,h' l n' 2 \C k \nin 2 ,hin 2 > = 
= V n k(n k - l)(5„^„ fc _ 2 + (2n k + l)(5 r ^„ fc + 

V {n k + l)(n k + 2)(5 n ^ lfc+2 . 

The calculation of the matrix element of H, the gen- 
erator of symmetry, gives us the following matrix: 

Z =< n^n!,, n^h'^H — H\ri\n 2 ,h\h 2 >= 
= uj{~(2cosh 2 9)5 n < ini 5 n2ri2 + 
+(ni + n 2 )S n ' xni S n2 n 2 }+ 

+ -^cosh 4 e{y/ni(ni - 1)<5„' 1 „ 1 _ 2 + (2ni + l)^/ ni + 

+ N /(m + l)(n 1 + 2)J„ ini+2 } (9) 

Diagonalizing the matrices R and Z numerically we ob- 
tain the energy eigenvalues of the Yang-Mills-Higgs sys- 
tem for the Hamiltonians H and H at finite-temperature. 
As it was mentioned the numerical energy levels de- 
pends on the dimension of the truncated matrix. Wc 
compute the numerical levels in double precision. The 
matrix dimension is 1156 x 1156, i.e. we calculate the 
first 1156 eigenvalues.. Then the statistical properties of 
the spectra are found. We use standard unfolding pro- 
cedure in order to remove the secular variation of the 
level density as a function of the energy E, for each value 
of the coupling constant the corresponding spectrum is 
mapped, by a numerical procedure described in |33| . 

One of the main characteristics of the statistical prop- 
erties of the spectra is the level spacing distribution 0, 
function. In this work we calculate the nearest-neighbor 
level-spacing distribution [2|, y, U, |7j . The nearest neigh- 
bor level spacings are defined as Si = E. L+ i — Ei, where 
Ei are the energies of the unfolded levels, which are ob- 
tained by the following way: The spectrum {Ei} is sepa- 
rated into smoothed average part and fluctuating parts. 
Then the number of the levels below E is counted and 
the following staircase function is defined: 

N(E)=N av (E)+N fluct (E). 
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The unfolded spectrum is finally obtained with the 
mapping 

Ei = Nav(Ei). 

Then the nearest level spacing distribution function 
P(S) is defined as the probability of S lying within the 
infinitesimal interval [S, S + dS] . 

For the quantum systems which are chaotic in the clas- 
sical limit this distribution function is the same as that of 
the random matrices 0)0- For systems which are regular 
in the classical limit its behaviour is close to a Poissonian 
distribution function. This distribution is usually taken 
to be a Gaussian with a parameter d: 

P(H) ~ exp{-Tr{HH + }/2d 2 ) : 

and the random matrix ensemble corresponding to this 
distribution is called the Gaussian ensemble. For Hamil- 
tonians invariant under rotational and time-reversal 
transformations the corresponding ensemble of matrices 
is called the Gaussian orthogonal ensemble (GOE). It 
was established @, LA L| that GOE describes the statis- 
tical fluctuation properties of a quantum system whose 
classical analog is completely chaotic. The nearest neigh- 
bor level spacing distribution for GOE is described by the 
Wigncr distribution: 

P(S) = InSexpi-^S 2 ). (10) 

The usual way to study the level statistics is to compare 
the calculated nearest-neighbor level-spacing distribution 
histogram with the Wigner distribution. 

For systems whose classical motion is neither regular 
nor fully chaotic (mixed dynamics) the level spacing dis- 
tribution will be intermediate between the Poisson and 
GOE limits. Several empirical functional forms for the 
distribution have been suggested for this case [4| . If the 
Hamiltonian is not time- reversal invariant, irrespective of 
its behavior under rotations, the Hamiltonian matrices 
are complex Hcrmitian and the corresponding ensemble 
is called Gaussian unitary ensemble (GUE). If the Hamil- 
tonian of the system is time-reversal invariant but not in- 
variant under rotations, then the corresponding ensemble 
is called the Gaussian symplectic ensemble. In the next 
section we will study the numerical results for the level 
spacing and then analyze them to classify them in one 
of the these categories. It is to be emphasized that our 
Yang-Mills-Higgs Hamiltonian is time-reversal and rota- 
tional invariant. 



IV. NUMERICAL RESULTS 

On diagonalizing the vacuum expectation value of the 
Hamiltonian H, physical observables, and H, the gener- 
ator of symmetry, given in Eq. 0and Eq. [51 respectively 
the spectra are analyzed by considering the level-spacing 
distribution. Then the criterion mentioned in the last 



section allows us to classify the system as chaotic or non- 
chaotic. In Figs. 1 and 2 we plot the level spacing distri- 
butions for different values of parameters u> and 9 for the 
energy spectrum calculated by diagonalizing of the ma- 
trix R. In Fig. 1 this distribution is given for the value 
u) = 0.92 for which the non-thermal case level spacing 
distribution is chaotic . As is seen from this figure 
for 9 = it is the same as the results of non-thermal cal- 
culations [23 . By increasing the temperature it becomes 
closer to a Gaussian distribution that means strength- 
ening of chaos in the thermal case. In Fig. 2 the level 
spacing distribution for u = 0.01 is plotted. For this 
value of oj the systems is regular at 6 = 0. However 
the increase of temperature leads to a chaotization of the 
system and P(S) becomes closer to the Gaussian distri- 
bution. Figs. 3 and 4 present level spacing distributions 
(for the same values of parameters as in Figs 1 and 2) for 
the spectrum calculated by diagonalizing of the matrix 
Z. By comparing Figs. 1, 2, 3 and 4, it is clear that 
at zero temperature level spacing distributions are the 
same for both methods. By increasing the temperature, 
the difference between spacing distributions for different 
methods becomes considerable. It is clear from Figs. 1 
and 2, that the level spacing distribution for ui = 0.92, 
9 = 0.01 is closer to the Gaussian distribution compared 
to ones for u> = 0.01, 9 — 0.01. The same behavior can 
be seen in histograms for u> — 0.92, 9 = 0.2 and uj = 0.01, 
9 = 0.2. The reason for such behavior can be understood 
from Fig. 5 where temperature T is plotted as a function 
of 9 for various values of to. It is clear that higher values 
of a; correspond to higher temperatures, while for smaller 
values of lu temperature is also small. 

However, for uj = 0.92 the difference between the re- 
sults for H and H is small while for uj = 0.01 there is 
considerable difference in the results even for 9 = 0.01 
(see Fig. 2 and 4). This indicates that the results from 
H, physical observables and H, generators of symmetry, 
are quite similar at high temperature while at low tem- 
perature the results are quite different. We know that at 
zero temperature H and H perform different functions. 
But the results indicate that at high temperature their 
chaotic behavior is similar. 

Thus in both cases increasing the temperature leads 
to a smooth transition from a Poissonian to a Gaussian 
form in the level spacing distribution. Furthermore, at 
higher temperatures both H and H lead to quite similar 
results. 



V. CONCLUSION 

Summarizing we have treated quantum chaos in gauge 
fields at finite temperature using a toy model, SU{2) 
Yang-Mills-Higgs system. To account for the finite - 
temperature effects we used the thermofield dynamics 
technique. The need for simultaneous exploration of the 
level fluctuations and the finite-temperature effects is dic- 
tated by recent advances in relativistic heavy ion collision 
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experiments, that allow the creation of hot and dense 
quark-gluon and hadronic matters |32j| . 

The lattice QCD calculations of hadronic matter and 
quark-gluon matter indicate that both systems exhibit 
strong chaotic dynamics 0, The present calcula- 

tions in a toy model appear to support such conclusions. 
Furthermore the study of the Yang-Mills-Higgs system at 
finite temperature establishes clearly that an increase of 
temperature of the system strengthens level fluctuations 
in the spectra. Then a transition from the Poissonian 
to a Gaussian level-spacing distribution does occur. It 
is to be anticipated that a study of the quark-gluon sys- 
tem in the relativistic heavy ion collisions at RHIC and 
later on at LHC would show the phenomenon of quantum 
chaos in a quantum field theoretic system. It is to be em- 
phasized that a proper study in (3 + 1) dimensions with 



full details of the Yang-Mills non-abelian field along with 
the Higgs scalar field is needed to make necessary con- 
clusions for a quark-gluon plasma. However the present 
study provides an indication of the possible outcome in 
more realistic studies. 
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